Semiclassical Accuracy in Phase Space for Regular and Chaotic Dynamics 
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A phase-space semiclassical approximation valid to 0{h) at short times is used to compare 
semiclassical accuracy for long-time and stationary observables in chaotic, stable, and mixed systems. 
Given the same level of semiclassical accuracy for the short time behavior, the squared semiclassical 
error in the chaotic system grows linearly in time, in contrast with quadratic growth in the classically 
stable system. In the chaotic system, the relative squared error at the Heisenberg time scales linearly 
with ftcff , allowing for unambiguous semiclassical determination of the eigenvalues and wave functions 
in the high-energy limit, while in the stable case the eigenvalue error always remains of the order of 
a mean level spacing. For a mixed classical phase space, eigenvalues associated with the chaotic sea 
can be semiclassically computed with greater accuracy than the ones associated with stable islands. 



I. INTRODUCTION 

Semiclassical methods have a long history traceable to 
the very beginnings of the "old quantum theory" and 
serve two interrelated purposes in many areas of physics. 
First, semiclassical methods provide valuable approxima- 
tion techniques in situations where a full quantum cal- 
culation is either impossible or unnecessary. Equally im- 
portantly, semiclassical methods provide a link between 
quantum results and our classical intuition, and allow us 
to separate physical behavior that is due simply to clas- 
sical paths and their interference from behavior that is 
attributable to non-classical processes such as tunneling 
or diffraction. 

For strongly chaotic systems, purely classical calcula- 
tions in d dimensions that ignore phase effects must break 
down at the mixing time or log time Tiog ^ A~^ IniV ^ 
X~^{d— 1) ln?i~g , where A is the maximal Lyapunov ex- 
ponent of the classical dynamics and N ~ ?i^ ~ ' is the 
effective dimension of the accessible Hilbert space, or the 
size of the accessible classical phase space in Planck cell 
units. This breakdown of c/assica^quantum correspon- 
dence occurs because beyond the mixing time, multiple 
classical paths connect a generic initial state to a generic 
final state, and interference effects become 0(1). On the 
other hand, in a series of papers. Heller and coworkers 
showed that semiclassical calculations in chaotic systems, 
which include the effect of interference between distinct 
classical paths, can follow the quantum propagator at 
times well beyond the mixing time [1]. An estimate for 
the breakdown time scale of semiclassical-quantum cor- 
respondence was obtained by quantifying the effects of 
caustics for a stadium billiard [2]. 

Over the past decade, significant light has been shed 
on the issue of semiclassical accuracy and its breakdown 
in diverse chaotic and regular systems. For example, 
Boasman has used a semiclassical approximation to the 
boundary integral method to obtain a semiclassical spec- 
trum for two-dimensional chaotic billiards, observing an 
overall semiclassical spectral shift as compared with the 
exact quantum spectrum, in addition to small random 
fluctuations [3]. On the other hand, Prosen and Robnik 
have shown the complete failure of torus quantization to 



reproduce the spectra of two-dimensional integrable bil- 
liards, such as the circle billiard [4], suggesting that inte- 
grability may in some cases lead to an increase of semi- 
classical errors; Rahav et al. have obtained more recent 
results consistent with this conclusion [5]. Primack and 
Smilansky were among the first to analyze semiclassical 
accuracy for three-dimensional chaotic systems, focusing 
on including corrections to the state-counting function 
beyond the leading Weyl term [6]. 

Main and collaborators have developed the power- 
ful harmonic inversion technique for accurate and effi- 
cient semiclassical calculations of energies, resonances, 
and matrix elements [7]. This technique, as well as the 
earlier cycle-expansion method [8] were applied recently 
to the four-sphere scattering problem, demonstrating a 
high degree of accuracy at a greatly reduced computa- 
tional cost compared with brute-force quantum calcula- 
tions [9]. Another promising recent approach, put for- 
ward by Vergini and coworkers, involves the accurate 
construction of quantum eigenstates as linear superpo- 
sitions of "scar functions" associated with short unstable 
periodic orbits [10]. 

In the time domain, statistical arguments concerning 
the propagation of semiclassical errors have shown that 
semiclassical error in chaotic systems accumulates inco- 
herently, and thus the squared error typically grows only 
linearly with time, in contrast with quadratic growth 
for the regular case [11]. Transforming to the the en- 
ergy domain, this implies that semiclassical methods 
are generically more accurate for computing wave func- 
tions and eigenvalues for chaotic systems than for reg- 
ular ones, in the S-eff — *■ (or high energy) limit. In 
particular, for d = 2, analytical arguments and numer- 
ical tests show that eigenvalues can be semiclassically 
resolved with great accuracy for chaotic systems, for suf- 
ficiently small Tics, while in the regular case even the or- 
der of eigenvalues cannot be unambiguously determined 
semiclassically. This result has been related to the reduc- 
tion of the quantization ambiguity in chaotic systems [12] 
and to the slower decay of fidelity in the presence of 
strong chaos (as long as the perturbation has non-zero 
diagonal matrix elements in the basis of the unperturbed 
system) [13]. 



Both the theoretical analysis and the numerical tests 
in Rcf. [11] were performed for scmiclassical evolu- 
tion in the position representation, i.e., for the Van- 
Vleck-Gutzwiller propagator [14]. Although well-suited 
for the model systems treated in that work, position- 
representation semiclassics suffers in general from the 
problem of proliferation of caustics, which eventually 
dominate the semiclassical propagator [2]. The prob- 
lem becomes particularly acute when one attempts to 
compare semiclassical dynamics in hard chaotic systems 
with that in a regular system or in a mixed phase space. 
Semiclassical calculations in a phase space basis are more 
natural from the point of view of classical-quantum cor- 
respondence and have the inherent advantage of allowing 
direct comparison between time evolution in chaotic, reg- 
ular and mixed systems, without the result being over- 
whelmed by the problem of position-space or momentum- 
space caustics. 

The aim of this paper is to improve our understanding 
of semiclassical accuracy in a phase space representation, 
as a function of time and hes, and to directly compare 
the behavior of the semiclassical error in chaotic, regular, 
and mixed systems. The organization is as follows. In 
Section II we briefly present the model and the method 
used for performing semiclassical and quantum calcula- 
tions in phase space. Theoretical expressions for semi- 
classical accuracy in chaotic systems are presented in Sec- 
tion III A, along with supporting numerical data from the 
model system. In addition to generalizing the analysis of 
Ref. [11] from position space to a phase space represe- 
nation, we exphctly test the prediction of Ref. [11] con- 
cerning the linear growth with time of the mean squared 
semiclassical error, as well as the prediction of linear de- 
crease with h of the error at the Heisenberg time. This 
is followed by a similar analysis for regular and mixed 
systems, in Sections IIIB and IIIC, respectively. We 
find that the behavior of the semiclassical error at the 
energy scale of a mean level spacing corresponds well to 
the effect of the quantization ambiguity on the spectrum, 
which was previously studied in Ref. [12]. Finally, Sec- 
tion IV summarizes the results and presents an outlook 
for the future. 



II. MODEL AND METHOD 

Although our theoretical analysis applies quite gen- 
erally to two-dimensional noninteracting systems, we 
simplify the numerical simulations by focusing on one- 
dimensional kicked maps, which share most scaling and 
other physical properties of this class of systems [15] . The 
discrete-time map can be regarded as a Poincare surface 
of section of a two-dimensional system with Hamiltonian 
dynamics. Specifically, we parametrize the one-step map 
on a toroidal classical phase space {q,p) G [0, 1) x [0, 1) 



where 



^o 



V{q) 



qi = qo +T'{pi) mod 1, 



(1) 






■ sin(27rg) 



1 K 

T{p) = -nipp^ + —^ sin(27rp) . 



(2) 



The dynamics is iterated to obtain classical evolution 
over many kicks (or many bounces in the correspond- 
ing two-dimensional Hamiltonian system). For values 
TOq — nip = 1 and < \Kq\, \Kp\ < 1, for example, we 
obtain the purely chaotic perturbed cat map, or kicked 



inverted oscillator, while for m„ 



-1, m,: 



1 and small 



Kg, Kp, the dynamics is predominantly regular, corre- 
sponding to a kicked regular oscillator. Parameters Kg 
and Kp are essential to introduce nonlinearity into the 
dynamics (if Kq ~ Kp = 0, the semiclassical propagator 
is exact, for any niq and nip). The above four parameters 
can also be adjusted to vary the Lyapunov exponent in 
the chaotic regime, or to study a mixed phase space, as 
we will see below in Section III. 

The one-step quantum evolution matrix for the above 
system takes the very simple form 



C/i = exp Hf{p)/h] exp HV{q)/h] , 



(3) 



Po 



Pi =Po - ^'(^o) mod 1 



which again may be iterated or diagonalized to obtain 
long-time or stationary behavior, Ut — [UiY- As dis- 
cussed in the introduction, we will apply this propagator 
to Gaussian wave packets (or coherent states) centered 
at phase space points (qk, Pk)'- 

(i>k{q) ^Uexp[-{q-qk)^/2h + ipk{q-qk)/fi], (4) 

where A/" is a normalization constant. 

Unfortunately, there is not a unique and universally 
used semiclassical approximation for wave packet evolu- 
tion, analogous to the Van-Vleck-Gutzwiller expression 
in position or momentum space. Several methods have 
been proposed that differ in both the order (in hes) of the 
semiclassical error at fixed time t and in the numerical 
size of that error. The so-called "thawed Gaussian" ap- 
proximation, for example, allows the shape of the Gaus- 
sian wave packet to change as it evolves under a locally 
quadratic potential V{q) [16]. An alternative approach 
uses "frozen" or unspreading wave packets [17]. Another 
coherent state method retains the stationary phase idea 
of the Van-Vleck-Gutzwiller propagator but extends dy- 
namics into complex phase space [18]. It is possible in- 
stead to work in complex time while retaining real initial 
conditions in phase space [19]. 

In the present work, we are not interested in reduc- 
ing the numerical size of the semiclassical error but only 
in understanding its scaling properties with t and hcs, 
for regular and chaotic systems. For this reason, we will 
choose what is a convenient method for our purposes, 
noting that the results would hold for any semiclassical 



approximation valid to the same order in hcs- We es- 
sentially use a variation of the thawed Gaussian method, 
extended to next-to-leading order in V^eff , and then cal- 
culate semiclassically the overlaps of the time-evolved 
"thawed" Gaussians with the Gaussians in the original 
basis [20]. 

Specifically, we start with a (non-orthogonal) complete 
set of A^ = l/2TTh Gaussians (f)j of the form given in 
Eq. (4), with the center points {qj,Pj) offset slightly from 
a rectangular grid to reduce numerical instabilities. The 
semiclassical overlap matrix 



Aa{j,k) = {(pj\(f>k] 



sc 



(5) 



is obtained analytically by Gaussian integration. To eval- 
uate the i— step semiclassical propagator At (j, k) between 



initial Gaussian (hi- and final Gaussian 



^j' 



we find real 



classical trajectories from {qo,Po) to {qt,Pt) in time t that 
minimize (go - Qk)'^ + (Po - Pk)'^ + [qt - qjY + {pt - PjY, 
i.e., all trajectories that start near the center of Gaus- 
sian k and end near the center of Gaussian j after t 
steps. Of course for fixed t and sufficiently small h 
(t < Tiog ^ A^^ \ogh~^), there will be at most one such 
trajectory, and in principle that is all we need even for 
our long-time analysis, as will be seen below. In practice 
however, for finite values of H we include all contribut- 
ing trajectories. For each trajectory, the potential V{q) 
is expanded to third order around the starting position 
of the trajectory, go- When this potential is applied to 
the original Gaussian (f>k , we obtain a wave packet of the 
form 

exp [a + b{q - qa) + c{q - qoY + d{q - q^f] (6) 
= exp [a + b(q - qo) + c{q - go)^] 
X [1 + d{q - qaf + 0{h)] , 

where a, 6, c, and d are complex numbers of order h~ , 
and therefore q — qo is 0{h^'^). We note that an "ex- 
tended" semiclassical dynamics [21], which truncates the 
expansion of the Hamiltonian at third order rather than 
second order is needed to keep the error in the one-step 
phase space propagator at 0{K), consistent with the er- 
ror in the Van-Vleck-Gutzwiller propagator in position 
space [14]. 

The wave packet of Eq. (6) may now be rewritten, via 
Fourier transform, as a momentum space wave packet 
having the same form but expanded in powers oi p ~ pi 
instead oi q — q^. The kinetic term T{p) of the Hamil- 
tonian may now be applied, again expanded to third or- 
der \n p — pi. Then, the packet is Fourier transformed 
back to position space and the procedure is repeated t 
times. At the end of t steps, we may analytically find the 
overlap between the semiclassically evolved t— step wave 
packet 4'k,sc(t), still having the form of Eq. (6), and the 
final Gaussian wave packet (f>j to obtain the semiclassical 
propagator At{j,k). If several classical paths lead from 
the vicinity of (j)k to the vicinity of (f>j in time t, their 



contributions must be summed to produce the semiclas- 
sical amplitude At{j, k), just as in the Gutzwiller expres- 
sion. As we will see in Section III A, for a chaotic system 
the long-time semiclassical propagator may be arbitrar- 
ily well approximated (in the h ^ limit) using only 
the matrix At for 1 <IC t <IC Tlog, where at most one path 
contributes to each matrix element. However, as we are 
dealing with finite h in our numerical simulations, we will 
always use the sum over all classical paths in numerical 
calculations. 



III. SEMICLASSICAL ACCURACY 
A. Chaotic Dynamics 

As discussed previously in the context of position- 
space semiclassical propagation, direct comparison be- 
tween quantum and semiclassical evolution at long times 
for a chaotic system, or between quantum and semi- 
classical stationary properties for such a system, faces 
the obstacle of the exponential proliferation of classical 
paths [1]; an analogous problem of exponential growth 
in the number of periodic orbits exists in the energy do- 
main [14,8,23]. This proliferation seemingly makes long- 
time semiclassical propagation in a classically chaotic 
system an exponentially harder problem than the full 
quantum evolution, puts into question the convergence 
of long-time semiclassical dynamics to any stationary be- 
havior, and prevents the comparison of semiclassical and 
quantum stationary properties for small h. The three- 
fold difficulty can be addressed using the idea that the 
Heisenberg uncertainty principle washes out information 
on scales below h, and thus the total amount of semiclas- 
sical information is finite for all times and scales only as 
a power of h. We can therefore collect, consolidate, and 
iterate semiclassical amplitude on sub-?i scales, obtaining 
the full semiclassical long-time dynamics to arbitrary ac- 
curacy in polynomial computation time. This "semiclas- 
sical path consolidation" idea has previously been used 
successfully to investigate long-time semiclassical accu- 
racy in the position representation for chaotic dynam- 
ics [20,11] and to demonstrate the semiclassical nature of 
dynamical localization in one dimension [22]. Conceptu- 
ally, the approach is similar to cycle expansion methods 
in periodic orbit theory [8,23]; however, no information 
about periodic orbits is needed here. Instead of account- 
ing for long-time semiclassical behavior in terms of peri- 
odic paths up to period Tpcriodic ~ 2]og, we use all short 
paths up to length t ^ 1. In the following, we adapt 
the methods of Ref. [11] to a phase space representation, 
and refer the reader to that earlier paper for a detailed 
discussion. 

We begin by noting that although semiclassical dynam- 
ics is not multiplicative, due to the fact that a concatena- 
tion of two stationary paths is in general not stationary, 
we may nevertheless write 



= [A,,AoiA,JO-,fc) + 0(?i), (7) 



where the 0(K) error is due to the intermediate sums 
being done exactly rather than by stationary phase, and 
the inverse of the semiclassical overlap matrix Aq is nec- 
essary due to non-orthogonality. In general, we may ap- 
proximate the true time-i semiclassical propagator At 
by evaluating the exact semiclassical dynamics to some 
"quantization time" r and then iterating the resulting 
matrix: 



At 



At mod T [^0 ^1 



1 \tlr\ 



(8) 



where [</r] is the integer part of t/r. We may call At.j 
the "r-semiclassical" propagator. For r = 1, At^T is the 
one-bounce semiclassical quantization pioneered by Bo- 
gomolny [24]. For a continuous-time system, the t <^\ 
limit is equivalent to quantum propagation via the Feyn- 
nian path integral approach. The exact time-i semiclas- 
sical propagator, on the other hand, is recovered in the 
opposite limit when the quantization time approaches t: 



A, 



t,T 



At as T ^ t . 



(9) 



In Ref. [11], it was shown analytically and numerically 
that the error \At^T{j,k) — Ai(j, /c)|^ falls off as Tc\/t in 
a chaotic system, where Tc\ is the time scale of classical 
correlations. This implies that for r ^ Td the approx- 
imate semiclassical correlator At,r is closer to the exact 
semiclassical correlator than either is to the quantum dy- 



\At,r - At\ <^ \At - Ut\ 



Thus 



\At 



utr^^^iAt 



Ut 



(10) 



allowing for an unambiguous determination of the error 
in the true semiclassical dynamics At at time t using At,r 
and permitting a study of the breakdown of the semiclas- 
sical approximation at long times t where performing an 
exact sum over O (e) classical paths is impractical or 
impossible. 

To confirm the convergence of the iterated propagator 
At,r to the true long-time semiclassical propagator At for 
semiclassical dynamics in phase space, and specifically 
the convergence of the semiclassical error in accordance 
with Eq. (10), we first compute, as a function of time t, 
the average r— semiclassical error defined as: 

Et^r = \\At^r - UtW" - Tr[At^r - Ut]HAt,r - Ut] 

^J2\^tAj,k)-Utij,k)\\ (11) 

The results are shown in Fig. 1 for a chaotic kicked map 
defined by parameters mq — nip — 1 and Kq — Kp = 



1/2, with semiclassical parameter N — 256. We notice 
the relatively poor agreement between the iterated semi- 
classical calculation for r = 1 and similar calculations 
for larger quantization times t. We also note that using 
the iterated propagator with short quantization time t 
overestimates the true size of the semiclassical error. At 
the same time, we observe rapid convergence of Et^r as 
T ^ 1, with the T = 5 and r = 6 curves lying almost on 
top of one another. Thus, the r-semiclassical error Et^r 
appears to be rapidly approaching the true semiclassical 



Et^J2\MJ,k)^Ut{j,k)\' 



(12) 



jife 




FIG. 1. The mean squared r-semiclassical error Et,T in a 
chaotic system is plotted as a function of t/N for semiclassical 
parameter N = 256 and for several values of the quantization 
time r. For r S> 1, Et,T is a reliable proxy for the true 
semiclassical error Et- The classical system parameters for 
Eq. (2) are m, — nip = 1; Kq — Kp — 1/2. The lower and 
upper solid curves represent the theoretical predictions for Et 
given by Eqs. (13) and (14), respectively, with Ci = 0.017 
and C2 = 0.037. 

We are now ready to investigate the semiclassical er- 
ror Et as a function of time t and semiclassical parameter 
N . For a chaotic system, we may assume that the errors 
associated with the semiclassical approximation add in- 
coherently as long as the times at which the errors occur 
are separated by at least the classical time scale Tci [H]- 
Since the squared error in the semiclassical approxima- 

we 



tion over one time step is Ei = 0{h?) 
obtain 

Et = Cih^t 



0{l/N' 



' N 



(13) 



where Ci ^ Td is a system-dependent constant and we 
take t — I to correspond to one period of the kicked 



map. The linear growth of the error predicted by Eq. (13) 
breaks down at times comparable to the Heisenberg time, 
where we must include an additional error term that is 
diagonal in the eigenbasis of the true quantum propaga- 
tor Ui [11]. The error associated with diagonal matrix 
elements adds coherently, leading to quadratic growth of 
the cumulative error in time. However, the fraction of 
diagonal matrix elements scales as /i = 1/A^. Eq. (13) 
must therefore be modified to read 






^h 



Ci 



(14) 



The data in Fig. 1 for r > 4 shows good agreement 
with the prediction of Eq. (14), which is indicated by 
the upper solid curve. The linear growth indicated by 
Eq. (13), shown as the lower solid line, is valid for times 
t short compared with the Heisenberg time TV. 



1 and 2. The scaling of the error at the Heisenberg time 
determines the feasibility of semiclassically computing in- 
dividual eigenstates and eigenvalues in the limit of small 
hcff, corresponding physically to the high-energy limit 
E ^ Egs- Based on Eq. (14), we predict the error at the 
Heisenberg time to be proportional to h: 



Et=N ^h[C\+ C2] 



(15) 



This prediction is tested in Figure 3, where the black 
squares represent the numerical data and the correspond- 
ing solid line is a best fit to a power-law form, Et=N — 
ah^ — aN~^, with /3 « 0.8. This is to be compared with 
the asymptotic prediction /3 = 1 for /i -^ 0. The falloff in 
the error with N shows that individual eigenstates and 
eigenvalues may be determined with ever improving accu- 
racy as iV — > 00. As we will find in the following section, 
this is in contrast with the situation for systems with 
regular classical dynamics (see also the white squares in 
Fig. 3). 



N=256 

N=144 

N=6 4 




chaotic dynamics 
regular dynamics 
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FIG. 2. The mean squared r-semiclassical error Et,T 
for a chaotic system is plotted as a function of t/N for 
r = 5 S> 1 and several values of the semiclassical parame- 
ter N = 1/h — 64, 128, 256. The classical system parameters 
are the same as in the previous figure. The lower and upper 
solid curves represent the theoretical predictions of Eqs. (13) 
and (14), respectively. The dotted line indicates the pre- 
dicted growth of the error for a system with regular dynamics, 
Et ^ t (see Section IIIB), and is shown to emphasize the 
qualitatively different behavior. 

In Fig. 2, we confirm the behavior predicted by 
Eqs. (13) and (14) as we vary the semiclassical parameter 
N = l/h. In this figure, the error Er,t has been scaled 
by a factor of N to make the curves at different values 
of N approximately coincide and to emphasize that the 
error at a fixed fraction of the Heisenberg time is falling 
off as /i ~ 1/A^ in the semiclassical limit h ^ 0. 

Specifically, we may ask about the size of the semiclas- 
sical error at the Heisenberg time itself, i.e. at t/N — 1, 
which corresponds to the right edge of the graph in Figs. 



FIG. 3. The semiclassical error at the Heisenberg time, 
Et=N,T is plotted for t = 5 ':^ 1, and for three values of 
the semiclassical parameter A'^ — l/h. Black squares cor- 
respond to the chaotic system of the previous two figures, 
while white squares correspond to the regular system of Sec- 
tion IIIB. The straight lines are fits to the power-law form 
Et=N ~ aN~^ , with the best fit giving /3chaotic = 0.8 and 
/Sroguiar = 0.1, to be Compared with the theoretical predic- 
tions /^chaotic = 1 (Eq. (15)) and /^regular = (Eq. (21)). 

The semiclassical spectrum and semiclassical eigen- 
states can be obtained in principle by computing the 
semiclassical propagator At for long times and transform- 
ing into the energy domain. However, since the semiclas- 
sical propagator At at long times becomes approximately 
multiplicative [11], 



A{t + l)KA^,A{t) 



(16) 



for some constant matrix y4*, it is much more conve- 
nient to diagonalize A* directly to obtain the semiclassi- 
cal eigenvalues and wave functions. We emphasize that 



A^ is neither the quantum evolution matrix Ui nor the 
semiclassical evolution matrix Ai for one time step, but 
is instead the effective one-step semiclassical propagator 
that describes semiclassical evolution at long times, and 
thus the stationary behavior of the semiclassical dynam- 
ics [11]. In practice, we may obtain A^, as the limit 



A^ = lim A,_7- 

T — *00 



lim A{t + 1) [A{t)]~^ 



(17) 



As discussed in Ref. [11], the convergence A^,^r -^ A^, 
is exponentially fast in r, at least for the position space 
semiclassical propagator: 



||^*,r A.^ 



,2„-At 



(18) 



In Fig. 4, we verify this convergence in the case of the 
phase space semiclassical propagator for two different val- 
ues of iV = l/h (white and black circles). The rate of 
convergence A is consistent with the classical value of the 
Lyapunov exponent, and is independent of h. The white 
triangles correspond to an example with a larger Lya- 
punov exponent {mq = 2, nip = 1, Kq ~ Kp ~ 1/2 in 
Eq. (2)), where the convergence with t is correspondingly 
faster. 
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FIG. 4. The convergence of the finite-time approximation 
to the long-time one-step semiclassical propagator A, is stud- 
ied for several systems and different values of the semiclassical 
parameter A'' = Xjh (see Eq. (17)). The circles represent data 
for the same system parameters that were used in the pre- 
vious three figures, white and black circles corresponding to 
A'' = 64 and N — 144, respectively. White triangles repre- 
sent data for A'^ = 64 with an alternative set of parameters: 
niq — 2, nip = 1, Kq — Kp = 1/2 in Eq. (2), having a larger 
Lyapunov exponent. The solid lines for the two systems are 
the predictions of Eq. (18), with A = cosh"^ (3/2) = 0.96 and 
A — cosh^^ (2) — 1.32, respectively. The white and black 
squares (A'^ = 64 and N = 256, respectively) represent data 
for the regular dynamics discussed in Section IIIB, where no 
convergence with r is predicted or observed. 



Exponentially fast convergence to A^ with r implies 
that the semiclassical spectrum and semiclassical wave 
functions can be obtained with very high accuracy using 
semiclassical dynamics for t > 1 but still short compared 
to the Heisenberg time t = N oi even the log time T^i^. 
In other words, all the information needed to calculate 
long-time or stationary semiclassical properties is already 
contained in the short-time classical behavior, well before 
before interference effects become relevant. 

The stationary semiclassical spectrum and wave func- 
tions can now be compared with their quantum ana- 
logues. From the linear scaling with h of the error in the 
time evolution at the Heisenberg time, Eq. (15), which 
has been tested above in Fig. 3, we can deduce that the 
mean squared error in the eigenvalues must also scale 
linearly with h, ignoring a possible overall shift in the 
spectrum [12] which is absent in the present system due 
to symmetry. Thus, 



F 



1 



N , 



A2 



1 

TV' 



(19) 



where the e^ and Cj^sc are the quantum and semiclas- 
sical eigenvalues, and A is the mean level spacing. In 
practice, this improvement in the semiclassical approxi- 
mation for individual eigenvalues as ft, ^ is difficult to 
measure due to numerical errors. For example, for the 
same chaotic system discussed previously {niq — mp — 1 , 
Kq^Kp = 1/2), F is already 1.3 • 10"^ for N = 36. 



B. Regular Dynamics 

We may easily change parameters in Eq. (2) to obtain 
fully or almost fully stable classical dynamics and then 
repeat the semiclassical calculations and analysis of Sec- 
tion III A. We choose nip = 1, niq = —1, Kp = Kq = 0.1. 
The small nonlinearity parameters Kp and Kq have been 
selected to reduce the semiclassical error in the short-time 
propagator; as we will see below, the semiclassical error 
grows much faster with time here than in the chaotic 
case. 

In a system with regular dynamics, a typical classi- 
cal trajectory repeatedly visits the same regions of phase 
space, and errors in the semiclassical approximation are 
expected to add coherently [12]. Thus, in contrast with 
the chaotic case, the squared difference between the time 
evolution matrix for quantum dynamics and its semiclas- 
sical counterpart is expected to grow quadratically with 
time: 






(20) 



where C is a classical constant that depends on the non- 
linearity of the system, as well as on the typical number 



of kicks needed for a typical classical trajectory to re- 
turn to the vicinity of its starting point. This quadratic 
growth of the error, even at times short compared to the 
Heisenberg time N, is to be contrasted with the result 
of Eq. (14) for a fully chaotic system. The prediction of 
Eq. (20) is tested in Fig. 5, where the quadratic growth is 
confirmed as well as the predicted scaling with the semi- 
classical parameter N = l/h. Furthermore, the growth 
of the semiclassical error with time is completely differ- 
ent in the regular and chaotic case, as can be seen from 
the dotted lines in Figs. 2 and 5. 
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FIG. 5. The mean squared r-semiclassical error Et,T for 
a system with regular dynamics is plotted as a function of 
t/N for r = 5 3> 1 and several values of the semiclassical 
parameter A*' = l/h = 64, 128, 256. The classical system 
parameters are nip — 1; ruq = — 1; Kq — Kp = 0.1. The solid 
curve represents the theoretical prediction of Eq. (20), with 
C = 0.0016. The dotted line corresponds to linear growth of 
the error with time, Et ^ t, applicable to the chaotic case only 
(see Eq. (13) and Fig. 2), and is shown to emphasize the very 
different scaling behavior in the case of regular dynamics. 



For a regular system at the Heisenberg time t — N , 
obtain an ?ioff-independent semiclassical error 



we 



for a two-dimensional system with regular classical dy- 
namics, the semiclassical error in the eigenvalues is pro- 
portional to the mean level spacing, implying that even 
the order of eigenvalues in the spectrum cannot be un- 
ambiguously determined using semiclassical methods. 

The problem, however, is more serious still, as the 
semiclassical dynamics for a regular system does not in 
general approach a stationary behavior at long times. 
We recall that for a chaotic system, the dynamics at 
long times approaches multiplication by a constant ma- 
trix A,: , whose eigenvalues and wave functions determine 
the stationary properties of the system. In contrast, for 
a regular system, the convergence of Eq. (17) docs not 
hold, since the Lyapunov exponent vanishes. This lack of 
convergence is observed in the squares plotted in Fig. 4, 
where it is seen that successive approximations to A^ 
differ from one another at 0{1/N^) — 0{h^ff). In other 
words, the eigenvalues of the matrix defining semiclassi- 
cal evolution from time t to t + 1 and the eigenvalues of 
the matrix defining semiclassical evolution from i -I- 1 to 
t + 2 differ from one another on the scale of a mean level 
spacing, so no unique semiclassical spectrum can be de- 
fined that describes the long-time semiclassical behavior. 

We note that a system with regular dynamics may be 
separable, in which case one may have a special set of co- 
ordinates for which semiclassical dynamics is exact (just 
as semiclassics may be exact for special chaotic systems 
such as the cat maps). The above results apply to the 
general situation where separability may not hold, e.g., a 
pseudo-integrable system or a generic polygonal billiard, 
and also to the separable case when the quantization is 
done in a set of coordinates other than the ones for which 
the equations of motion separate. Assuming the semiclas- 
sics is not exact, and independent of the initial size of the 
semiclassical error, the semiclassical accuracy will pro- 
gressively improve in the hcs ^ or high-energy limit as 
long as the Lyapunov exponent A is nonzero, until even- 
tually individual eigenvalues and wave functions become 
semiclassically resolvable. In the case of zero Lyapunov 
exponent, this improvement does not occur. 



Et=N — C . 



(21) 



to be contrasted with the diminishing semiclassical error 
at the Heisenberg time in the h^s ^^ limit for a chaotic 
system, as indicated by Eq. (15). The Heisenberg-timc 
error for our regular system is plotted for several values of 
N in Fig. 3. We note that the Heisenberg-time semiclas- 
sical error is larger for the regular system as compared 
with a chaotic system at the same value of N , despite 
the fact that the one-step semiclassical error is larger in 
the chaotic case. 

The 0{h ) error in the semiclassical evolution at the 
Heisenberg time, as indicated by Eq. (21), suggests that 
semiclassical eigenvalues and wave functions, if they ex- 
ist, do not approach the corresponding quantum eigen- 
values and wave functions in the ficS — *■ limit. Instead, 



C. Mixed Dynamics 

Generic two-dimensional systems are neither fully reg- 
ular nor fully chaotic, and it is therefore of interest 
to study the issue of semiclassical-quantum correspon- 
dence in the general regime of "soft chaos." A mixed 
classical phase space can be obtained using parameters 
TUq ~ Kp = 0, TOp = Kq = 1 in Eq. (2); for this system 
approximately 48% of phase space is associated with the 
chaotic sea and the remainder consists of stable islands. 
Based on our discussion in Sections III A and IIIB on 
the very different behavior of semiclassical accuracy in 
chaotic and regular systems, respectively, it is natural to 
ask whether semiclassical accuracy may vary with initial 
conditions in the case of a mixed phase space. 



We define a local version of the mean squared eigen- 
value error introduced in Eq. (19): 
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(b) 



where (j)k is one of the Gaussian wave packets introduced 
in Section II, ipi and Cj are the eigenstates and eigenval- 
ues of the quantum dynamics, and ei,sc are the semi- 
classically obtained counterparts to e^. In other words, 
F^i^ measures the error in the semiclassical eigenvalues, 
weighing each eigenvalue error by the overlap of the cor- 
responding eigenstate with 0^ . A contour plot of F^,, ver- 
sus phase space coordinates qk, Pk is shown in Fig. 6 (a), 
for N =l/h = 256. 



We see that the semiclassical error is peaked in the 
major stable regions of phase space, particularly in the 
large stable island surrounding the q = p = stable fixed 
point, and to a somewhat lesser extent in the islands as- 
sociated with the period-2 orbit a,t p = 1/2. In contrast, 
F^^, remains low in the region of the chaotic sea, for ex- 
ample, in the vicinity of the unstable orbit at q — 1/2, 
p = 0. The contour plot in Fig. 6 (b) shows the fraction 
of each wave packet 4>k consisting of stable trajectories, 
and the similarity between the main features in the two 
parts of the figure strongly suggests a correspondence 
between semiclassical accuracy and classical phase space 
structure. 




(a) 



FIG. 6. (a) The weighted semiclassical eigenvalue error F^,. 
is plotted as a function of phase space location {qk,Pk) for 
semiclassical parameter N = 256. The contour curves corre- 
spond to F = 0.001, 0.003, 0.005, 0.007 (the thickest curve 
indicates the largest error). The semiclassical eigenvalues are 
obtained by diagonalizing At,o in Eq. (17). (b) For each wave 
packet (pk used in (a), the fraction of stable trajectories for 
that wave packet is calculated classically and again plotted 
as a function of wave packet location. The contour curves 
correspond to stable fractions of 0.6, 0.9, 0.975 (the thickest 
curve corresponding to the most stable region). 

The total semiclassical error for a mixed system is of 
course dominated by the error associated with the stable 
regions, and scales in the same way as the error for a 
regular system in Section III B. 



IV. SUMMARY 

Phase-space semiclassical propagation allows us to 
make direct comparison of semiclassical validity in 
chaotic and stable classical systems. Using the same 
semiclassical approximation in both cases results in a 
semiclassical error that scales with Ti in the same way 
at short times. However, the growth of the error with 
time is very different in the two situations. In the regular 
case, the error grows coherently because each trajectory 
repeatedly visits the same regions of space phase; the 
mean squared error therefore grows quadratically with 
time. In the chaotic case, this coherence effect does not 
occur at times short compared with the Heisenberg time, 
resulting in a linear growth of the mean squared error. 

At the Heisenberg time itself, the mean squared error 
in the propagator matrix elements becomes 0(1) in the 
case of a classically stable dynamics, making it impossi- 
ble in general to speak of well-defined semiclassical wave 
functions or eigenvalues, i.e. ones that are independent 
of the choice of semiclassical coordinates. For a given 
choice of coordinates, semiclassical quantization generi- 



cally will produce eigenvalues differing by the order of 
a mean level spacing from their quantum counterparts. 
Different semiclassical quantizations of a regular system 
will produce spectra differing from each other at the same 
order, making it impossible to uniquely determine even 
the order of eigenvalues in the spectrum via semiclassical 
methods (unless a particularly favorable set of coordi- 
nates can be chosen where semiclassics happens to be 
exact, e.g., for a separable dynamics). 
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In contrast, semiclassical dynamics at the Heisenberg 
time for a classically chaotic system becomes increasingly 
accurate as the system energy is increased. In the energy 
domain, the semiclassical error becomes a progressively 
smaller fraction of a mean level spacing, so the spectrum 
can be semiclassically determined with arbitrarily high 
accuracy when very highly excited states are considered. 
The convergence of semiclassical to quantum behavior 
for chaotic system is expected to be independent of the 
particular semiclassical method chosen (for example, it is 
independent of whether a position, momentum, or phase 
space semiclassics is used) as long as the methods have 
the same scaling with h at fixed time. 

All calculations in the present paper were performed 
for time-dependent one-dimensional maps, whose scal- 
ing properties are equivalent to those of two-dimensional 
Hamiltonian systems. In d — 3 dimensions, or in an in- 
teracting system, the Heisenberg time grows as a higher 
power of h" than in the two-dimensional single-particle 
case, resulting in a larger accumulated semiclassical error 
by the Heisenberg time for both chaotic and regular sys- 
tems. For example, the same scaling argument that leads 
to Eq. (15) for d = 2 chaotic systems predicts 0(1) semi- 
classical errors at the Heisenberg time for chaotic sys- 
tems, independent of energy, i.e. eigenvalue errors that 
remain a constant fraction of a mean level spacing. In 
other words, the breakdown time of the semiclassical ap- 
proximation will be proportional to the Heisenberg time 
in three dimensions, even when the dynamics is chaotic 
(and much shorter than the Heisenberg time for regular 
dynamics) . 

For d > 4, e.g., in the case of two interacting par- 
ticles in two dimensions with no conserved quantities 
apart from total energy, the semiclassical approximation 
is expected to break down well before the Heisenberg 
time, even when the dynamics is fully chaotic. It would 
be interesting to investigate this behavior quantitatively 
for model systems, and also to ascertain how a higher- 
order semiclassical approximation may enable semiclas- 
sical methods to remain valid for interacting systems. 



[1] S. Tomsovic and E. J Heller, Phys. Rev. Lett. 67, 664 
(1991); M. A. Sepulveda, S. Tomsovic and E. J Heller, 
Phys. Rev. Lett. 69, 402 (1992). 
[2] S. Tomsovic and E. J Heller, Phys. Rev. E47, 282 (1993). 
[3] P. A. Boasman, Nonlinearity 7, 485 (1994). 
[4] T. Prosen and M. Robnik, J. Phys. A 26, L37 (1993). 
[5] S. Rahav, O. Agam, and S. Fishman, J. Phys. A 32, 7093 

(1999). 
[6] H. Primack and U. Smilansky, Phys. Rep. 327, 1 (2000); 
J. Phys. A 31, 6253 (1998); see also T. Prosen, Phys. 
Lett. A 233, 323 (1997). 
[7] J. Main, Phys. Rep. 316, 233 (1999); K. Weibert, J. 

Main, and G. Wunner, Phys. Lett. A, 289, 329 (2001). 
[8] P. Cvitanovic and B. Eckhardt, Phys. Rev. Lett. 63, 823 

(1989). 
[9] J. Main, G. Wunner, E. Atilgan, H. S. Taylor, and P. A. 
Dando, Phys. Lett. A 205, 176 (2002). 
[10] E. G. Vergini, J. Phys. A 33, 4709 (2000); E. G. Vergini 

and G. G. Carlo, tbtd. 33, 4717 (2000). 
[11] L. Kaplan, Phys. Rev. E bS, 2983 (1998). 
[12] L. Kaplan, New J. Phys. 4, 90 (2002). 
[13] T. Prosen and M. Znidaric, J. Phys. A 35, 1455 (2002); 
T. Prosen and T. H. Seligman, ibid. 35, 4707 (2002); T. 
Prosen and M. Znidaric, New J. Phys. 5, 109 (2003). 
[14] M. Gutzwiller, J. Math. Phys. 12, 343 (1971). 
[15] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. 
Rev. Lett. 49, 509 (1984); A. Altland and M. R. Zirn- 
bauer, Phys. Rev. Lett. 77, 4536 (1996). 
[16] E. J. Heller, J. Chem. Phys. 62, 1544 (1975). 
[17] M. F. Herman and E. Kluk, Chem. Phys. 91, 27 (1984). 
[18] R. G. Littlejohn, Phys. Rep. 138, 193 (1986); D. Huber, 
E. J. Heller, and R. G. Littlejohn, J. Chem. Phys. 89, 
2003 (1988). 
[19] W. H. Miller and T. F. George, J. Chem. Phys. 56, 5668 

(1972). 
[20] F.-M. Dittes, E. Doron, and U. Smilansky, Phys. Rev. E 
49, R963 (1994); L. Kaplan and E. J. Holler, Phys. Rev. 
Lett. 76, 1453 (1996). 
[21] A. K. Pattanayak and W. C. Schieve, Phys. Rev. E 50, 

3601 (1994). 
[22] L. Kaplan, Phys. Rev. Lett. 81, 3371-4 (1998). 
[23] G. Tanner, P. Scherer, E. B. Bogomolny, B. Eckhardt, 

and D. Wintgen, Phys. Rev. Lett. 67, 2410 (1991). 
[24] E. B. Bogomolny, Comm. At. Mol. Phys. 25, 67 (1990). 



